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Abstract: Before effective vaccines become available, antiviral drugs are considered as the 
major control strategies for a pandemic influenza. However, perhaps such control strategies 
can be severely hindered by the low-efficacy of antiviral drugs. For this reason, using 
antiviral drugs and an isolation strategy is included in our study. A compartmental model 
that allows for imported exposed individuals and asymptomatic cases is used to evaluate 
the effectiveness of control strategies via antiviral prophylaxis and isolation. Simulations 
show that isolation strategy plays a prominent role in containing transmission when antiviral 
drugs are not effective enough. Moreover, relatively few infected individuals need to be 
isolated per day. Because the accurate calculations of the needed numbers of antiviral 
drugs and the isolated infected are not easily available, we give two simple expressions 
approximating these numbers. We also derive an estimation for the total cost of these 
intervention strategies. These estimations obtained by a simple method provide a useful 
reference for the management department about the epidemic preparedness plans. 
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1. Introduction 

The influenza virus, such as avian influenza H5N1 virus, is a great threat to public health [1-3]. 
Because we do not know what new influenza strain will appear in the future, it is possibly more 
dangerous and more difficult to contain its quick spread than any known influenza strain [4]. Vaccines 
are recognized as an effective means to prevent the spread of pandemic influenza; however, it takes 
about 180 days to develop an effective vaccine when a new strain of pandemic influenza arises [3,5,6]. 
Before effective vaccines become available, antiviral drugs, such as Tamiflu, and measures that reduce 
exposure to infected cases are used to mitigate the transmission of influenza virus [3,7]. To be specific, 
these emergency measures consist of avoiding close contact as much as possible, isolating infectious 
individuals in the hospital or at home, wearing personal protective equipment, such as masks and gloves, 
and the use of antiviral drugs for prevention and treatment [4,7]. 

Some previous work discussed the use of antiviral drugs for the control of influenza. 
Longini et al. [8] explored the effectiveness of targeted antiviral prophylaxis, in combination with other 
control measures. Black et al. [9] studied household-based interventions with the use of the antiviral 
drugs for the control of pandemic influenza. Antiviral drugs have been proven to be effective against 
currently circulating strains of pandemic influenza, so many countries set up a stockpile of antiviral 
drugs for preparedness plans [10-14]. However, the efficacy of antiviral drug against a new pandemic 
influenza strain is not yet known, and it is difficult to confirm whether current antiviral drug is effective or 
to what extent it is effective, against a novel influenza virus strain [14] . If the reproductive number cannot 
fall below one by using antiviral drugs, there will be a widespread outbreak of the disease. Therefore, it 
is very necessary to implement other control measures in conjunction with the use of antiviral drugs. 

Several previous studies have shown that the isolation strategy is an effective control measure. Sunmi 
Lee et al. [15] used optimal control theory to demonstrate that the isolation strategy has a prominent 
role when antiviral resources are limited. Yan et al. [16] showed that using isolation and quarantining 
strategies as much as possible at the start of the outbreak are very critical. 

Based on these considerations above, in this manuscript, we have incorporated the use of 
antiviral drugs and the isolation strategy in our model. We use a simple compartmental model to evaluate 
the effect of using antiviral drugs for prophylaxis and the isolation strategy to reduce transmission. It is 
impractical that antiviral drugs would be used only for prophylaxis and not for the treatment of cases, so 
we take into account antiviral prophylaxis supposing that antiviral drugs are used for treating index cases 
first, as Merler et al. [17]. We give two simple expressions, which approximate the needed numbers 
of antiviral drugs and the isolated infected, respectively, providing a practical reference for the sizes 
of antiviral stockpiles and healthcare resources (such as beds and personnel) for the possible epidemic 
in advance. Specifically, we consider the cost of these intervention strategies before the time when an 
effective vaccine against the new strain becomes available. We also give a concise expression, which 
approximates the total intervention cost. 
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2. Methods 

2.1. The Basic Model 

In this section, we introduce the influenza transmission model studied in this paper. This model 
is based upon a SLIAR model, which is an extended version of the standard SEIR model and is 
developed by Arino et al. [18], but includes the importing cases from other regions and interventions 
with isolation and the use of antiviral drugs. Figure 1 summarizes the model schematically. We 
consider a homogeneous population whose individuals mix uniformly and suppose that the importation 
approximately equals the exportation, neglecting births and natural deaths, so that the population remains 
almost unchanged. For simplicity, the mortality associated with the disease has not been included in our 
model explicitly, and the qualitative behavior of the model is not changed by this simplification [3]. 

Figure 1. The summary of an influenza transmission model, where A = f m (elQ +ijv) +crA, 
(3 is the transmission rate. Exposed individuals leave the compartment at the rate k. A 
proportion 1 — p of infected individuals becomes asymptomatic cases and recovers in I/73 
days, and the rest of infected individuals develop symptoms. A fraction u of symptomatic 
infected individuals is isolated and recovers in I/71 days; the remainder of symptomatic 
cases recover in 1/72 days. The rates of importing newly exposed persons and asymptomatic 
cases from other regions are a u and a 2 t, respectively. The transmission rate can be reduced 
from /3 to f3f m because of the use of antiviral drugs. 




We write S for the percentage of the people who are susceptible, L for the fraction of exposed 
individuals (i.e., latent individuals — infected, but not infectious), Iq for the proportion who is 
symptomatic cases and has been isolated, I N for the proportion who is symptomatic cases, but not 
isolated, A for the fraction of the population who are asymptomatic, but infectious and R for the 
proportion who has recovered. The transmission and recovery are modeled by the following equations: 

S = —/3AS 

L = (3 AS -kL + a u 

Iq = upkL - 7 1 J Q 

I N =p(l - u)kL - j 2 In 

A= (l-p)kL-~f 3 A + a 2 t 

R = 71 I Q + 72 I N + J3A 
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with A = / m (eJg + In) + crA (where f m , e and a are interpreted below) and initial conditions: 

S(0) = S o , L(0) = L 0 , I Q (0) = I QO , I N (0) = I m , A(0) = A 0 , R(0) = R 0 

Susceptible individuals proceed to the latent compartment at the rate {3(f m (elQ + In) + aA)S, where 
(3 is the transmission rate and S denotes dS/dt. We assume that infectivity and symptoms of influenza 
start at the same time [19], and each infected person who has obvious symptoms presents to the health 
service when his or her symptoms appear. As Becker and Wang [14], we suppose that each newly 
confirmed case causes the distributing of m doses of antiviral drugs to the case and persons who have 
close contact with that case. These drugs are used for treating the confirmed index case and providing 
prophylaxis to the associates of this index case [8,17]. For simplicity, in the following, the words antiviral 
prophylaxis mean that antiviral drugs are used for treating the confirmed index case and providing 
prophylaxis to the contacts of this index case. Since distributing the first few doses of antiviral drugs to 
the closest associates of the index case is probably more efficient in reducing transmission, as Becker and 
Wang [14], we adopt the following form of f m , which means that the transmission rate can be reduced 
from (3 to (3f m by dispensing m antiviral doses per case, where: 

f m = a + (l-a)exp(-bm) (2) 

for various values of a and b satisfying 0 < a < 1 and b > 0. Obviously, f m decreases from 1 to a with 
m increasing. 

Exposed members leave the compartment at the rate k, and a proportion p of exposed members 
develops symptoms, whereas the remainder proceed to an asymptomatic infectious compartment. A 
proportion u of symptomatic cases have been isolated and go to the recovered class at the rate 7 1? while 
the rest of symptomatic cases recovered at the rate 72. The infectivity of isolated cases is reduced by a 
factor of e, with 0 < e < 1. The case e = 0 corresponds to the most perfect scenario when isolated cases 
have almost no contact at all with the susceptible individuals and the transmission rate is reduced to 
zero. The case e > 0 represents the situation when contacts between isolated and susceptible individuals 
are not avoided. Gumel et al. [20] interpreted the modification parameter e as the level of hygiene 
precautions during isolation. In fact, parameter e tunes the effectiveness of the isolation strategy, and 
it can somehow mimic the effect of a delay in the implementation of the interventions. The infectivity 
of asymptomatic individuals is reduced by a factor of a, with 0 < a < 1 . In addition, asymptomatic 
individuals progress to the recovered compartment at the rate 73. For simplicity, we assume that the 
transmission rate (3 and the recovery rates 71, 72, 73 are not affected by seasonal variation. 

In addition, we suppose that symptomatic cases will not travel and that latent individuals are 
potentially traveling, because they have no symptoms [21]. We also assume that asymptomatic cases 
are potentially traveling. The model (1) includes an inflow of exposed individuals into the community at 
a rate a lt and an inflow of asymptomatic cases into the community at a rate a 2 t- Although the rate 
of importing exposed individuals and the rate of importing asymptomatic cases are time-dependent 
in practice, in this manuscript, we will concentrate on constant importation rates, i.e., a u = «i 
and «2t = a 2 - 
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2.2. Reproductive Numbers 

In epidemiological models, the basic reproduction number is one of the most important quantities, 
denoted by R 0 [22,23]. It is the average number of infections produced by an infectious individual 
in a wholly susceptible population when no public health interventions are implemented [23,24]. In 
the presence of control measures, it is called the control reproduction number [23]. In our model, 
the control measures include antiviral prophylaxis and isolation, and we denote the corresponding control 
reproduction number by R c . Applying the approach developed by van den Driessche and Watmough [25] 
to our model equations shows that the next-generation matrix K is: 



K 



7i 72 73 

0 



e/3fm 




(3a 


7i 


72 


73 


0 


0 


0 


0 


0 


0 


0 


0 


0 



(3) 



0 

y 0 0 0 0 ) 

Since the rank of this matrix is 1, its spectral radius and trace are equal, which gives: 

flt = Wm (!»! + £&Z!0) + f(izrt ] and R^^L^JLlA) (4 ) 

71 72 73 72 73 

2.3. Intervention Cost for an Outbreak 

Interventions can contain the spread of an outbreak, but they are costly. One of our aims is to derive 
an approximate expression for the intervention cost for an epidemic outbreak. Then, firstly, we need to 
model the cost of interventions. 

In this paper, we concentrate on two major intervention measures: antiviral prophylaxis and isolation. 
Following Ludkovski and Niemi [26], the total intervention cost is equal to a fixed cost plus a variable 
cost. A fixed cost is paid once when the intervention begins, and a variable cost is paid each time period 
and may depend on the latest number of actual infected individuals and isolated infected individuals. 
For simplicity, we assume the variable cost of antiviral prophylaxis is in proportion to the cumulative 
number of individuals who have received antiviral drugs up to time t v , i.e., 

Ci / mpkL(t)dt (5) 
Jo 

where t v is the time when an effective vaccine against the novel strain of influenza can be supplied, the 
coefficient c x is the cost per dose of antiviral drug and the integral expression: 

'* £ v 

pkL(t)dt (6) 



o 

represents the cumulative percentage of individuals who have received antiviral drugs up to time t v . The 
fixed cost of antiviral prophylaxis is denoted by cs(t v ), which is a one-time upfront cost. Similarly, the 
variable cost of the isolation strategy is in proportion to the cumulative number of the infected who have 
been isolated up to time t v , i.e., 

c 2 / upkL(t)dt (7) 



o 
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where the coefficient c 2 is the cost of one infected individual who has been isolated during his or her 
infection period and the integral expression: 

ty 

upkL(t)dt (8) 

o 

represents the cumulative percentage of the infected individuals who have been isolated up to time t v . 
The fixed cost of the isolation strategy is denoted by c±(t v ), which is a one-time upfront cost. 
Then, the total intervention cost for an outbreak on [0, t v \ is: 

pt v pt v 
ci I mpkL(t)dt + c 2 / upkL(t)dt + c 3 (t v ) + c^{t v ) (9) 



JO JO 

For simplicity, we take c 3 (t v ) and c 4 (t t) ) to be constants and denote c 3 (t„) = c 3 , c±(t v ) = c 4 

2.4. Approximating the Needed Number of Antiviral Doses Provided to Community Members, the 
Number of Isolated Infected Individuals and the Intervention Cost 

Obviously, by numerical calculations, we can obtain the needed number of antiviral doses provided 
to community members, the number of isolated infected individuals and the intervention cost from 
Equations (1), (6), (8) and (9). It cannot be denied that the computational process is somewhat complex, 
and it is not intuitive and convenient for the management department. It is more instructive and 
convenient to work with concise expressions that approximate these numbers. This is because we can 
know clearly what factors influence these numbers and how these factors influence them, through these 
simple expressions. Next, we consider the approximate values of these numbers under the assumption 
that the epidemic can be contained. 

When R c is always smaller than 1, the depletion of susceptible individuals can almost be ignored, in 
other words, S ~ 1. Inspired by the idea from [14], we may rewrite: 

L w f3[f m (eI Q + I N ) + a A] - kL + a x (10) 

By solving the simultaneous equations formed by the Equation (10) and the third to fifth equations 
in the model (1) with letting the right side of each equation be zero, we obtain the equilibrium 
value of the percentage of exposed individuals. The percentage of exposed individuals is close 
to its equilibrium value of («i + (3 a 0:2/73) /( k (1 — R c ))- Replacing L{t) with its equilibrium 
value («! + /3cra 2 /73)/(fc(l — R c )) in Equation (6), we can obtain that the approximate cumulative 
percentage of individuals who have received antiviral drugs up to time t v is pt v (ai + (5o a 2 /73)/(l — R c ) ■ 
Therefore, the number of antiviral doses that need to be distributed to community members on [0, t v ] 
is approximately: 

mpt v Ba 

=-(«! + — a 2 ) (11) 

1 - Rc 73 

and the cost of antiviral prophylaxis is approximately: 



(12) 
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Similarly, by using Equation (8), we can obtain with the same methodology that the percentage of 
cumulative isolated infected individuals is approximately: 

upt v f Ba 

a.\ H a 2 ) (13) 



1-Rc 73 
and the cost of the isolation intervention is approximately: 

uptv ( da 

c 2 - ^"(«iH «2) + c 4 (14) 

1 - Re 73 

Then, the total intervention cost for an outbreak on [0, t v ] is approximately: 

~(+ \ m P tv ( I Ba \ upt v Ba 

c{t v ;m,u) = ci- — (ai H a 2 ) + c 2 - — («i H a 2 ) + c 3 + c 4 (15) 

1 - -Re 73 1 - it c 73 

2.5. Parameters 

Table 1 shows baseline values of all parameters of the model, where the values of a±, a 2 , 8, 71, 72, 
73 and k are rates per day. We assume that a baseline initial value of R 0 = 1.5, as in [14]. The latency 
period for influenza is supposed to be equal to the incubation period and is set at 1.5 days on average, 
namely k = 1/1.5 [17]. We also suppose that the possibility that an infected individual is developing 
symptoms of influenza is 0.5, namely p = 0.5 [17]. We further assume that symptomatic cases are twice 
infectious as those without symptoms, namely a = 0.5 [8,17]. The recovery rate 7; = 1/1.5 (i = 1, 2, 3) 
means that the length of the infectious period is assumed to be 1.5 days [17,21]. The parameter e = 1/6 
assumes that the infectivity of isolated cases drops to 1/6 of the infectivity of symptomatic cases who are 
not isolated. The value t v = 180 days supposes that an effective vaccine will take 180 days to develop, 
manufacture and dispense, and a = 0.6 assumes that the effect of antiviral drugs is to reduce the risk of 
transmission per close contact by no more than 40%. In addition, we suppose that each person is initially 
susceptible, i.e., Sq = 1 and L 0 = Iqo = I no = A 0 = R 0 = 0. The epidemic is seeded by the imported 
infectious individuals. Unfortunately, there are no good data on the cost coefficients c\, c 2 , c 3 and c 4 ; 
therefore, we give roughly the relative values of the cost coefficients. 

Table 1. Baseline values for model parameters. 
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3. Results 

3.1. Antiviral Drugs Prophylaxis and Isolation Strategy 

Because effective vaccines will take about six months to produce once a novel influenza virus has 
been confirmed, the use of antiviral drugs is one of the most important intervention measures in the case 
of a pandemic [4,23]. We first consider that control strategy only via antiviral drugs for prophylaxis 
during an outbreak of pandemic influenza. 

Figure 2. (a) The plot of the control reproduction number R c as a function of doses m, 
where u = 0 and other parameters assume the values of Table 1 ; (b) the plot for which R c is 
equal to one. Every parameter point (m, u) that corresponds to R c < 1 lies above this curve 
and R c > 1 for those points that lie below the R c = 1 curve. 



(a) (b) 




Figure 2a shows the plot of the control reproduction number R c as a function of doses m, where 
u = 0 and other parameters suppose the values of Table 1 . We see that the control reproduction number 
R c remains greater than one, even for substantial values of doses m. If R c is above one, the disease 
will widely spread in the population; hence, the successful containment of the epidemic is to reduce the 
reproduction number below one. It is impossible to contain the epidemic only by the use of the antiviral 
drugs if antiviral drugs have low-efficacy in reducing infectiousness against a new virus strain. Therefore, 
it is very necessary to implement other control measures in conjunction with the use of antiviral drugs. 

Isolation is one of the most effective methods to contain the spread of an outbreak [27]. For our 
purposes, we consider that the control strategy is taking the isolation strategy in conjunction with the use 
of antiviral drugs for prophylaxis. There are two possibilities when taking antiviral prophylaxis and the 
isolation strategy: R 0 > R c > 1 and R 0 > 1 > R c . The first corresponds to failure containment, and 
infected arrivals could lead to an outbreak. In this paper, we focus on the second scenario, i.e., R c can be 
reduced to below one when taking antiviral prophylaxis and the isolation strategy widely and effectively. 

Figure 2b shows the values of m and u for which the control reproduction number R c is equal to one, 
namely R c = 1 for every parameter point (m, u) that lies on the curve. Meanwhile, every parameter 
point (m, u) that corresponds to R c < 1 lies above this curve and R c > 1 for those points that lie below 
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this curve. By comparing the two curves in Figure 2, we can see that the reproduction number that is 
greater than one can be reduced to below one by the implementation of isolation of cases. Suppose the 
stockpile of antiviral drugs is sufficient; at least 20% of the infected individuals need to be isolated for 
reducing the reproduction number below one. On the other hand, at least 60% of the infected individuals 
need to be isolated if no antiviral drugs are provided. 

To provide a more intuitive result, Figure 3 shows plots of the cumulative percentage of infected 
individuals up to time t and the percentage of isolated infected individuals versus time t, respectively, as 
predicted by the model (1) with the parameter values of the model listed in Table 1, for different propor- 
tions of u, where u = 0, 0.3, 0.5, 0.7, 0.9 and m = 15. Comparing the curves in Figure 3a, we can see that 
taking antiviral prophylaxis combined with the isolation strategy can reduce substantially the cumulative 
number of infected individuals. At 180 days, the cumulative number of infected individuals when 
u = 0.3 drops to 1/25 of the cumulative number of infected individuals when no cases are isolated. By 
way of illustration, for a city with one million population members, allowing for one imported infective 
individual and one imported exposed individual per day, the cumulative number of infected individuals 
will be less than 2, 134 and the number of daily isolated cases will be less than six for a duration of six 
months if half of the symptomatic cases can be isolated. 

Figure 3. (a) Plots of the cumulative percentage of infected individuals up to time t; and 
(b) plots of the percentage of isolated cases versus time t, as predicted by the model (1) with 
parameters: R 0 = 1.5, a = 0.6, b = 0.2, m = 15 and other parameters supposing the values 
of Table 1, for different values of u, where u = 0, 0.3, 0.5, 0.7, 0.9. 
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In addition, from Figure 3a, we can see that the cumulative percentage of infected individuals 
is always a monotonically decreasing function of u, i.e., increasing the proportion of isolation can 
reduce the number of infected individuals; but the percentage of isolated infected individuals is only a 
monotonically decreasing function of u, after some fixed time point. This can be interpreted as follows. 
With the doses m fixed, R c decreases when the proportion u of isolation increases. Hence, in a short 
time horizon, the percentage of the isolated cases increases as u increases; but in a long time horizon, as 
R c decreases, the percentage of cumulative infected individuals decreases, and so does the percentage of 
the isolated cases. In conclusion, for containing the epidemic better, we need to keep R c under a level as 
low as possible, which implies that it is important to isolate infected individuals as much as possible in 
the early stage of the disease. 

From Figure 3b, we can see that as long as a few infected individuals need to be isolated every day 
and the reproductive number is reduced to below one, then the total number of the infected individuals 
can be reduced greatly (see Figure 3a). 

It is obvious from the above that implementing an isolation strategy has an obvious effect on the 
containment of the epidemic, especially if antiviral drugs have low-efficacy in reducing infectiousness. 

3.2. Comparisons between the Approximate Values and the Actual Values of the Percentage of 
Cumulative Symptomatic Cases, the Percentage of Cumulative Isolated Infected Individuals and the 
Intervention Cost 

Figure 4a,b shows the actual values and the approximate values of the cumulative percentages of 
symptomatic cases and isolated cases as functions of time t when R c = 0.8. Comparing the dashed 
lines with the solid lines in Figure 4, we can find that the dashed lines are located above the solid lines 
all of the time, but there is little difference between them, and they have the same change trends. This 
illustrates that the difference between the actual values and the approximate values is very little. More 
specifically, for a duration of 180 days, for every million individuals, the absolute error of the cumulative 
number of symptomatic individuals is 129, and the absolute error of the cumulative number of isolated 
cases is 80. Thus, in practice, we can take the approximate values as the estimations of the actual 
values, and the errors are quite small. On the other hand, the expressions of estimated values are simple 
and easy to be calculated. Based on the above merits, the approximate method is very useful for the 
management department to make control strategies. However, it should be noted that the absolute errors 
of the cumulative percentages of symptomatic individuals and isolated cases increase with the control 
reproduction number R c increasing. In other words, the smaller the control reproduction number R c is, 
the more accurate the approximation is. As we can see in Figure 4c, the curve begins relatively flat and 
after a while becomes very steep when R c exceeds 0.9. This implies that the approximation precision 
is relatively high if the reproduction number can be reduced below 0.9, and the approximation precision 
is significantly lowered when R c is over 0.9. The low accuracy may be caused by the approximation of 
S ~ 1. The deviation between the percentage of susceptible individuals and one cannot be ignored when 
R c is close to one. 
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Figure 4. (a,b) Plots of the actual values and the approximate values of the cumulative 
percentages of infected individuals and isolated cases as functions of time t, respectively, 
where R c = 0.8, m = 15 and u = 0.62; (c) plots of the absolute errors of the cumulative 
percentage of symptomatic cases as a function of R c ; (d) plots of the absolute errors of the 
cumulative percentage of isolated cases as a function of R c . 




R R 

c c 

Figure 5 shows the actual intervention cost and the approximate intervention cost as given by the 
model (1) with c\ = l,c 2 = 50, c 3 = 100, c 4 = 100 and other parameters supposing the values of 
Table 1. The three graphs all show that the difference between the actual intervention cost and the 
approximate intervention cost is very little. The reason is that the approximate values of the cumulative 
percentage of symptomatic cases and the cumulative percentage of isolated cases are very close to their 
actual values, and the cost is a linear combination of them, so the result is inevitable. 

At the same time, it is of interest to note that in Figure 5b, the actual intervention cost is not necessarily 
a monotonically increasing function of doses m, and there is a turning point in the actual intervention 
cost curve, where it is minimum. A similar situation takes place in the approximate intervention cost 
curve. For all parameter values, the actual intervention cost and the approximate intervention cost 
show a trend downward at first, when they reach their bottoms, and then begin to rise again, as m 
increases. This phenomenon tells us that in order to minimize the cost of intervention measures, a modest 
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number of doses is enough, and too lager values of doses of antiviral drugs to reduce transmission only 
waste doses. 



Figure 5. (a) Plots of the actual intervention cost and the approximate intervention cost as 
functions of time t, where R c = 0.8, and the cost coefficients c\ = 1, c 2 = 50, c 3 = 100 
and C4 = 100; (b) plots of the actual intervention cost and the approximate intervention cost 
as functions of doses m, where R c = 0.8 and the parameters as (a); (c) plots of the actual 
intervention cost and the approximate intervention cost as functions of it, where m — 15 and 
the parameters as (a). 



200.1 



200.05 



200 



200.1 



200.05 



200 



200.08 

CO 

8 200.06 

CD 

jE 200.04 
200.02 



(a) 



■ The actual cost 



The approximate cost 




80 100 
Time (days) 

(b) 



180 



1 The actual cost 
The approximate cost 




10 



15 



20 



25 



30 



35 



40 



45 



50 



m 

(c) 



■ The actual cost 



The approximate cost 



0.4 



0.45 



0.5 



0.55 



0.6 



0.65 0.7 
U 



0.75 



0.8 



0.85 



0.9 



0.95 



These results have obvious implications for a public health service system. First, we give an esti- 
mation of the number of the cumulative isolated infected, so that the hospital can prepare the resources 
(such as beds and personnel) for the possible epidemic in advance. Second, the estimation of the required 
number of antiviral drugs provides a reference for the management department about the antiviral stock- 
pile size for controlling the epidemic. This is very significant for the management department, because 
a small stockpile of antiviral drugs cannot satisfy the demand of the containment epidemic, whereas, 
a very large stockpile of antiviral drugs is a great waste of medical resources. Finally, the estimation 
of the cost of intervention will help the government with financial preparation prior to the possible 
pandemic influenza. 
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3.3. Sensitivity Analysis 

The baseline values of the parameters in Table 1 may not be suitable to a newly emerged strain of 
influenza virus. Therefore, we carry out a sensitivity analysis to investigate the effect of varying these 
parameters, such as the basic reproduction number (Ro), the infectivity reduction factor for the isolated 
individuals (parameter e), the proportion of symptomatic cases (parameter p), the relative infectivity of 
asymptomatic cases (parameter a), the imported rate of exposed individuals (parameter a{), the imported 
rate of asymptomatic cases (parameter a 2 ), the efficacy of antiviral drugs (parameters a and b), the 
infectious period (parameters I/71, I/72, I/73) and the coefficients of the intervention cost (parameters 
ci, C2, C3 and C4). Figures 6-13 show their impact on the results in Sections 3.1 and 3.2. When conducting 
the sensitivity analysis for a parameter, we assume the values of other parameters as Table 1, unless 
indicated otherwise. 

Figure 6. (a) Plots of the cumulative percentages of infected individuals and isolated cases 
during 180 days, as functions of R 0 , where m — 15, u — 0.9; (b) plots for which R c is equal 
to one for different values of R 0 , where R 0 = 1.4, 1.5, 1.7, 2.0. Every parameter point (m, u) 
that corresponds to R c < 1 lies above this curve and R c > 1 for those points that lie below 
the R c = 1 curve. 



(a) (b) Rc=1 




The effect of the basic reproduction number: Earlier in the manuscript, we supposed that the basic 
reproduction number R 0 is 1.5. Figure 6a shows the effect of changing this assumption on the cumulative 
percentage of infected individuals. The cumulative percentage of infected individuals is quite sensitive 
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to this change. For fixed parameters m and u, the cumulative percentage of infected individuals increases 
as Rq increases, and it increases rapidly when R 0 surpasses some value. This is because the reproduction 
number cannot be reduced to below one for those values of R 0 that exceed some value. The change of R 0 
gives the same impact on the cumulative percentage of isolated cases. Figure 6b shows plots for which 
R c is equal to one for different values of R 0 , where R 0 = 1.4, 1.5, 1.7, 2.0. Every parameter point (m, u) 
that satisfies R c < 1 lies above the corresponding curve. The increase in the reproduction number R 0 
makes R c = 1 curve upward and shrinks significantly the set of parameter points for which R c < 1. This 
implies that it needs more effort to contain the transmission with larger _R 0 . Moreover, the containment 
of the epidemic will not succeed if the reproduction number is too high, such as R 0 > 2.5. 

Figure 7. The effect of the infectivity reduction factor for isolated individuals, (a) Plots of 
the cumulative percentages of infected individuals and isolated cases as functions of e, for 
a duration of six months, where m = 15, u = 0.9; (b) plots for which R c = 1, for various 
values of e, where e = 1/6, 1/3, 1/2, 2/3. 
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The effect of the infectivity reduction factor for the isolated individual: The above results assume 
that the infectivity reduction factor (parameter e) for the isolated individuals equals 1/6. The impact of 
altering this assumption is shown in Figure 7. The change of e has a significant impact on the cumulative 
percentages of infected individuals and isolated cases (see Figure 7a). The low value of e results in a 
much less cumulative percentages of infected individuals and isolated cases, while larger e generates 
much more infected individuals and isolated cases. Especially, the slopes of the two curves suddenly 
grow larger when the value of e exceeds 0.8, which means that the effectiveness of isolation strategies 
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decreases significantly for e > 0.8. This implies that the effectiveness of the isolation strategies is 
very sensitive to the parameter e. Since the earlier the isolation begins, the better the effectiveness of 
the control strategies is, isolation strategies should begin as early as possible for better effectiveness of 
interventions. To obtain successful containment of the epidemic, the objective of the implementation of 
control measures is to bring R c below one. Figure 7b shows the curves for which R c is equal to one for 
different values of e = 1/6, 1/3, 1/2, 2/3. As the parameter e increases, the set of scenarios for which 
containment is achievable becomes smaller and smaller. Even containment will fail if the parameter e is 
larger than some value. 

The effect of the proportion of symptomatic cases: Figure 8 shows the effect of varying the proportion 
of symptomatic cases (parameter p). The effectiveness of control measures is very sensitive to the 
proportion of symptomatic cases. For the fixed parameters m and u, the cumulative percentage of 
infected individuals decreases substantially with the proportion of symptomatic cases increasing. This is 
because more infected individuals are treated/isolated and more susceptible individuals receive antiviral 
prophylaxis. In Figure 8b, the set of scenarios for which containment is successful shrinks significantly 
as the proportion of symptomatic cases declines. If the proportion of symptomatic cases falls below a 
certain value (for example p < 0.2), the epidemic is hard to control. 

Figure 8. The effect of the proportion of symptomatic cases, (a) Plots of the cumulative 
percentages of infected individuals and isolated cases as functions of p, for a duration of six 
months, where m — 15, u — 0.9; (b) plots for which R c = 1, for various values of parameter 
p, where p = 0.25, 0.35, 0.5, 0.75. 
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The effect of the relative infectivity of asymptomatic cases: The relative infectivity of asymptomatic 
cases (parameter a) is studied in the range 0.003-0.8. The impact of varying the parameter o on the 
cumulative percentages of infected individuals and isolated cases is shown in Figure 9a. Unsurprisingly, 
the effectiveness of control measures declines as the infectivity of asymptomatic cases increases. A low 
infectivity of asymptomatic cases (a = 0.003) corresponds to a smaller cumulative number of infected 
individuals; conversely, a high infectivity (such as a = 0.8) generates more infected cases. However on 
the whole, the effect caused by the change of a is not very significant, since the two curves in Figure 
9a always gently rise. Figure 9b shows the plots for which R c = 1, for different values of parameter a, 
where a = 0.003, 0.3, 0.5, 0.8. The set of scenarios for which containment is available becomes smaller 
as the values of parameter a grows larger, and the reproduction number can be reduced below one for all 
the values of a ranging from zero to one. 

Figure 9. The effect of the relative infectivity of asymptomatic cases, (a) Plots of 
the cumulative percentages of infected individuals and isolated cases as functions of the 
parameter a, for a duration of six months, where u = 0.9, m = 15; (b) plots for which 
R c = 1, for various values of parameter a, where a = 0.003, 0.3, 0.5, 0.8. 
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The effect of the imported rates of exposed individuals and asymptomatic cases: We use three values 
for the imported rates of exposed individuals and asymptomatic cases («j = 10~ 7 , 10~ 6 , 10~ 5 , i = 1,2) 
to illustrate the effects of varying the assumed values of them. The change of the imported rate of exposed 
individuals makes no difference to the control reproduction number R c , according to the expression 
of R c (Formula (4)). However, this change has an impact on the cumulative percentages of infected 
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individuals and isolated cases. It should be noted that for a fixed value of a 2 , the cumulative percentages 
of infected individuals and isolated cases increase linearly with the increase in the imported rate of 
exposed individuals (ati). From Figure 10b, we can see that the effect of varying the imported rate of 
asymptomatic cases (a 2 ) is similar to that of varying the parameter a\. 

Figure 10. The effect of the imported rates of exposed individuals and asymptomatic cases, 
(a) Plots of the cumulative percentages of infected individuals and isolated cases as functions 
of «i, for a duration of six months, where m — 15, u = 0.9; (b) plots of the cumulative 
percentages of infected individuals and isolated cases as functions of a 2 , for a duration of 
six months, where m — 15, u — 0.9. 
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The effect of the efficacy of antiviral drugs: In Section 2.1, we assume that the effectiveness of 
dispensing m antiviral doses per case is to reduce the transmission rate from j3 to f3f m , where 
f m = a + (1 — a) ex-p(-bm). The sensitivity of b has been discussed in the literature [14]. For large 
doses m, the effectiveness of antiviral drugs in reducing transmission largely depends on the parameter 
a. Therefore, we only consider the effect of varying the parameter a. The impact of the change of 
a on the cumulative percentages of infected individuals and isolated cases is shown in Figure 1 1 . The 
cumulative percentages of infected individuals and isolated cases increase as the effectiveness of antiviral 
drugs declines for the fixed values of m and u. It is interesting to note that the influence of the change 
of a on the cumulative percentages of infected individuals and isolated cases becomes smaller with the 
proportion of isolation increasing. Figure lid shows the curves for which R c = 1, for various values 
of parameter a = 0.3, 0.5, 0.6, 0.8. The set of scenarios for which containment is successful becomes 
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large with the parameter a decreasing. In other words, containment is more likely to be successful when 
the value of a is small. Unfortunately, there is no good data to estimate the antiviral efficacy for a new 
pandemic strain before a control policy is instituted. In this case, regardless of whether antiviral drugs 
can effectively mitigate the transmission of a new pandemic strain, non-pharmaceutical measures, such 
as case isolation, combined with antiviral drugs, can be used to reduce the spread of influenza strain. 

Figure 11. The effect of varying the efficacy of antiviral drugs, (a) Plots of the cumulative 
percentages of infected individuals and isolated cases as functions of the parameter a, where 
m = 15 and u = 0.5; (b) same as (a), but for m = 15 and u = 0.7; (c) same as (a), but for 
m = 15 and u = 0.9; (d) plots for which R c = 1, for various values of parameter a, where 
a = 0.3,0.5,0.6,0.8. 
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The effect of the infectious period: The above results are based on the assumption that the 
infectious period is 1.5 days [17,21]. However, some models suppose that the infectious period is four 
days [3,14]. We carry out sensitivity analysis on the value of the infectious period. The impact of varying 
infectious periods on the effectiveness of control strategies is shown in Figure 12. The cumulative 
percentage of infected individuals is a little sensitive to the change of the infectious period with e = 1/6, 
but it is very sensitive to the change of the infectious period when e = 2/3 (see Figure 12a). The 
cumulative percentage of isolated individuals has the similar result as the cumulative percentage of 
infected individuals. Figure 12c,d shows the curves for which R c = 1 for different infectious periods 
l/7i, [i = 1, 2, 3) with e = 1/6 and e = 2/3, respectively, where = 1.5, 2.5, 3.5,4.5(i = 1, 2, 3). 
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The set of scenarios for which containment is successful expands with the infectious period extending, 
and the expansion is more obvious for a larger value of parameter e. 

Figure 12. The effect of varying parameters determining infectious period, (a) Plots of 
the cumulative percentage of infected individuals as functions of the infectious period, for 
e = l/6,l/3,2/3, respectively ; (b) same as (a) , but for the cumulative percentage of isolated 
individuals; (c) plots for which R c = 1 for different infectious period 1/7$, (z = 1,2, 3) with 
e = 1/6, where 1/7; = 1.5, 2.5, 3.5, 4.5; (d) same as (c), but with e = 2/3. 
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The effect of the coefficients of the intervention cost: In order to make a sensitivity analysis on the 
coefficients of the intervention cost (parameters ci, c 2 ) by comparing the actual cost with the approximate 
cost, we vary the coefficients c x and c 2 , respectively. Figure 13 shows the actual intervention cost 
and the approximate intervention cost as functions of time t using four different coefficients (c± = 1, 
c 2 = 50; ci = l,c 2 = 200; c\ = 50, c 2 = 1; c\ = 200, c 2 = 1, respectively). The two cost curves in 
every graph are quite close to each other. This illustrates that the changes of c\ and c 2 could not have a 
great influence on the difference between the actual intervention cost and the approximate intervention 
cost. Obviously, the changes of c 3 and c 4 have little impact on the conclusion that there is little difference 
between the actual intervention cost and the approximate intervention cost. 
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Figure 13. Plots of the actual intervention cost and the approximate intervention cost versus 
time t for various values of c\ and c 2 , where c\ = 1,50, 200, c 2 = 1,50, 200, R c = 0.8, 
c 3 = 100 and c 4 = 100. 
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4. Discussion 

One of our objectives was to see whether taking the isolation strategy should be included in an attempt 
to contain a newly emerged pandemic when antiviral drugs have low efficacy in reducing infectiousness 
against a new virus strain. The conclusion is a strong recommendation that implementing an isolation 
strategy is very effective if antiviral drugs are not effective enough. 

At the beginning of an outbreak, the efficacy of antiviral drug against a newly emerged pandemic 
influenza strain is not yet known, and it is difficult to confirm whether current antiviral drug is effective 
or to what extent it is effective. In this case, the implementation of the isolation strategy gives a chance 
to restrict the impact of pandemic influenza by lowering infection rates. Therefore, the isolation strategy 
is always effective for the uncertainty of the efficacy of the antiviral drugs against a newly emerged 
influenza strain. In addition, sensitivity analysis implies that the effectiveness of the isolation strategies 
is very sensitive to the parameter e. Since the earlier the isolation begins, the better the effectiveness 
of the control strategies is, isolation strategies should begin as early as possible for better effectiveness 
of interventions. 

The large fraction of asymptomatic individuals is one of the key features of influenza; 
Hayward et al. [28] pointed out that no symptoms were present in three quarters of the infected 
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individuals for the 2009 H1N1 pandemic. We therefore consider those infected individuals who have no 
symptoms. The conclusion is that containment is not successful if the proportion of asymptomatic cases 
exceeds a specific value, even though the intervention strategies are efficient enough. In other words, 
if taking antiviral prophylaxis and the isolation strategy cannot contain transmission at the beginning 
of an outbreak, the most likely reason for this failure is that a large number of infected individuals 
fail to display symptoms. In this case, the most effective measure for containing the epidemic may 
be to isolate the symptomatic cases who present to the health services and place the closest associates 
of the symptomatic cases into self-isolation, because they are possible asymptomatic cases. In short, 
the government should focus on antiviral prophylaxis and the isolation strategy in the early stage of an 
outbreak, and after a while, if containment is not achievable and the strain continues to spread, other 
intervention measures (such as self-isolation [27]) should be considered. This problem will be further 
explored in future work. 

Because the accurate calculations of the needed numbers of antiviral drugs and the isolated infected 
are not easily available, we give two simple expressions approximating these numbers. We also derive an 
estimation for the total cost of these intervention strategies. These estimations that are gotten by a simple 
method provide a useful reference for the management department about the epidemic preparedness 
plans. These conclusions depend largely on the condition R c < 1 (which means successful containment), 
and these results can be applied to a wide range of disease characteristics. This methodology still can be 
used to approximate other intervention costs, only if the intervention strategy can reduce the reproduction 
number to less than one. 

In addition, in order to confirm the effectiveness of antiviral drugs against the novel influenza virus 
strain, we must collect some data at the beginning of the local outbreak. According to the result in [14], 
the data on 100-200 household outbreaks are needed to examine an antiviral effect on transmission. 

On the basis of a simple compartmental model, some practical results have been obtained. 
These results are very meaningful and important for policy makers in making intervention 
strategies. Unfortunately, the results are not applicable to an epidemic with a high proportion 
of asymptomatic cases. 
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